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Abstract 


An attempt is made to account for the #2 disturbances in middle latitudes 
associated with the geomagnetic storms in term of the vertical drift of the electron. 
For this purpose, at first it is shown that in middle latitudes the variations of foF2 
on the storm day are classified in two types, negative and positive disturbances, 
the former of which is the same as the representative variation in high latitudes 
while the latter the analogous to the variation in low latitudes (including the 
equator), on the other hand, the variations of h’F2 are analogous to those in high 
and low latitudes, regardless of the seasons. 

The vertical drift of the electron is considered to be caused by an electric 
field deduced from the disturbance-daily variation of the geomagnetic field. The 
results of the calculations for individual states show that two types of disturbances 
in middle latitudes seem to be explained as an effect of the vertical drift of the 
electron, though the coincidence between the observed and calculated for the posi- 
tive type is not sufficient. From the above results, together with the results 
previously obtained for the equatorial zone, it is found that the occurrence of the 
positive or negative disturbances is ascribed to the conditions that the phase of the 
drift velocity of the electron on the quiet day is the same as that of the drift velocity 
on the disturbed day (negative) or the phase of the former differs from that of the 
latter by about 180° (positive). It is also found that the seasonal and latitudinal 
variations of the #2 disturbance are due to the existence of the Sg dynamo current 


in the ionosphere. 
1. Introduction 


The ionospheric F2 disturbances associated with geomagnetic storms have been 
studied theoretically by several authors. The theories are divided generally into three 
kinds. The first of them is the thermal expansion theory [1], [2], the second is the 
electron drift theory [3], [4] and the third is the recombination theory [5]. Of these 
theories, the electron drift theory is considered to be the most suitable one, because 
the ionospheric disturbances have the following characteristics : 

(a) The disturbances in the F2 regin are caused almost simultaneously in many 
stations of the world. me 3 . 

(b) When the disturbance begins in the /2 region, disturbances begin simultane- 
ously in the Fl and even in the £ regions [6]. The daily variations and the seasonal or 
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latitudinal changes of the deviations from the normal of the critical frequencies of the 
F1 and E regions are significantly analogous to those in the £2 region. [7]. 

Based on the above characteristics, the explanation of the F2 disturbances in 
middle latitudes is at first attempted for an average state of the disturbances by the 
drift theory by the author [8]. The results of the study show that the effects of the 
electron drift are very important. However, it seems to be better to study the F2 dis- 
turbances for their individual state, because they are fairly different from the averge 
state. From this point of view, the drift theory is also apllied by the author to the 
study of the individual F2 disturbances in the equatorial zone, using the data at 
Huancayo. The results reveal that the individual #2 disturbances in the equatorial 
zone are well explained as the effect of the vertical drift of the electron [9]. 

In this paper it is attempted to explain the individual F2 disturbances in middle 
latitudes by the electron drift theory and furthermore, based on the results of the cal- 
culations for the F2 disturbances in middle latitudes and in the equatorial zone, the 
features of the variations of the F2 disturbances with latitude, season, and sunspot or 
geomagnetic activities are physically explained. For these purposes individual F2 distur- 
bances are classified into two kinds and the features of their seasonal or latitudinal 
occurrences are stitistically clarified. 


2. Some Statistical Results of /'2 disturbances 


Statistical studies of the F2 disturbances are carried out by many authors [1], [3], 
[10]-[16] and the storm time deviations of the foF2 and h’F2 (critical frequency and 
virtual height of the F2 region) and their variations with season, latitude, and sunspot 
or geomagnetic activities are discussed. However, as the F2 disturbances in middle 
latitudes are complicated compared with those in other latitudes, the results derived 
somtimes seem to differ depending on the period, the year of the data and the number 
of the geomagnetic storms used there. This is because in many studies the attention 
is not paid for the type of the F2 disturbances. For this reason it is better to classify 
the daily variation of the deriation of the F2 disturbance by a suitable method and 
take statistics of the respective types of the disturbances. Up to the present time, the 
classification of the F2 disturbances has been attempted by some authors [10], [11], but 
theoretical bases of them are ambiguous. In this section we attempt a classification 
of the F2 disturbances and make their features clear. 

Glancing over many disturbances over the world, it is found that the daily varia- 
tion of fF 2 is classified into the following two types: 

(1) The f\F2 decreases from the normal in daytime with a maximum depression 
about at noon and recovers to the normal value or rather increases at night. 

(2) The foF2 tends to increase in daytime or/and at night. Sometimes foF2 
decreases early in the morning. 

In an actual case, one of the above variations is superposed on the other. We 
call tentatively the type (1) the ‘negative disturbance’and the type (2) the ‘positive dis- 
turbance’ (compare with Appleton and Piggot [11]). These two types of the F2.disturb- 
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ances are basic types, independently of the height variation (it is known in the follwing 
that the height deviations are always the same in every latitudes) and the physical 
meaning of the occurence of them-is shown in the section 7. The mean daily variations 
of the deviations (4MF2) in two:types are shown in Fig. 1 for five stations, College 
(geographic latit. 64.°9, geomagnetic latit. 64.°5) Washington (39.°0, 50.°0), Watheroo 
(-30.°19,-41.°8), Kokubunji (35.°7, 26.°0), and Huancayo (-12.°03,-0.°6). The blanks in 
Fig. 1 show rarity of the occurrence of the positive disturbance in high latitudes and 
negative disturbance in low latitudes or in the equator. 
Next we treat the 
W’F2 variation on the College 
storm day. We _ should 


pay an attention to the Ls ello 
derivation of the deviation p VR, aar 
of h’F2, because of the in- / : 

clusion of the retardation a ate Se Washington 


effect of the radio wave 
through the Fl region, 
Especially, in the case of 
the negative disturbance Teme aw! yenenie 
this effect is very great. Bare ee Wotheroo, 


Positi Oa 


The evidences are shown 
in Fig. 2 and A and Din- 
Fig. 6. Since the foF2 gb. at 7 Wire 
becomes very low due to MM / oe 
the disturbance, the differ- ad 
ence between f,F2 and 
foF 1 (critical frequeney of | 
the Fl region) becomes  g,+ aii 
very small. Then the \¢ 
deviation of h’F2 becoms 
large abruptly, sometimes 
it reaches to infinity. 
When the f,F2 is far 
greater than the ffl the 
deviation of the height is 4 

602 4 6 8 10 i2 4 16 18 202 6 
very small and rather the ocak Tie Local Time 


height decreases from the Fig. 1 Observed mean daily variations of AfoF2 dotted line and 
Ah’ F2 (deviative values from the normal of foF2 and h’F2) in the 

normal value. The retar- negative and positive disturbances for five stations. The variation 
dation effect is absent at is the mean of the remarkable deviations during the years Oct. 
5 saddent 1956-Jan. 1957 for College, 1949-1956 for Washington, 1938-1944 
night, sets up suddenly for Wetheroo, 1949-1956 for Kokubunji and 1938-1944 for Huancayo. 


at sunrise and disappears 
gradually as the sunset approaches. This effect is included not only in the A’F2, 


=] 


4 


but also in the hApF2 (height of 
the maximum electron density of 
the F2 region), because the latter is 
deduced from the calculation as to 
the one parabolic region and the 
effect of the other region is not 
taken into account. 

We must sometimes pay an 
attention to the retardation effect in 
the positive disturbanc in a meaning 
somwhat different from the above. 
The reason is as follows. The retada- 
tion effect is included in the height 


Fe nn 
0) 24 hours 48 


Storm time 


Fig. 2 An example of the /2 negative disturbance. 
It begins at 19h, March 2, 1956. The time is the 
nearest local time before the beginning of the geo- 
manetic storm, and this time is taken as zero of the 


on a quiet day in summer in middle 
or higher latitudes, because the differ- 
ence between fF'2 and fF'1 is small. 
So that when the positive disturbance occurs 
in summer the retardation effect is reduced, 
therefore the deviation of the height from the 
normal value is negative. 
According to the theoretical investigations 
[17], [18], [19], retardation depends on the 
values of foF1/fF2 and fy/fF2, where fF2 is 
the frequency of the wave reflected from the 
F2 rgion and fg the gyro-frequency, and the 
larger these ratios, the larger is the retardation. 
In Fig. 3 the values of the deviations of h’F2 
(4h’F2) in many F2 disturbances recorded at 
Washington, Watheroo and Kokubunji are 
plotted against the value of f\F1/fF2. The 
line in the figure is the value of the retarda- 
tion of the height against of (F'1/fF2, deduced 
from the theoretical calculation, assuming that 
fin /fF2=1/4, that the distance between the 
height of the maximum electron density of the 
F1 region and the lowest height of the F2 region 
is 150km and that the electron density above 
the maximum density of the F1 region is the 
same as the maximum density of the F1 region. 
Comparing the theoretical values with the ob- 
served ones, the coincidence is good for the large 
values of foF1/foF2 and it is not good for the 


storm time. 


0 a2 0.4 08 


0.6 
KFI/tF2 
fF FF2 
Fig. 3 The magnitude of 4h’ F2 against the 
value of foll/foF2. The full line in the 
figure represents the calculated value against 
the value of foF1l/fF2 where fF2 is the 
frequency of the radio wave by which the 
virtual height is measured. 
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smaller values of it. This is because that (F2~/fF2 for fF1/f\r2~1 and fol ceria 
for the lower value of AF 1/f,F2, if in the present case we assume that the f /2 repre- 
sents the frequency by which the-virtual height is measured. 

In Fig. 1 the mean daily variations of the deviative values (4h’F 2) of the height are 
shown for five stations, after a correction for the retardation is made using the theore- 
tical value shown in Fig. 3. Correction is made mainly for the negative disturbance. 
As it is difficult to correct the data at College for the lack of /,.F'1, the height deviations 
during Oct. 1956-Jan. 1957, whence the difference between foF1 and fyF2 is considered 
to be large, are taken. From Fig. 1 it is found that the daily variation of the deviative 
value of the height is almost the same for the stations from equator to high latitudes, 
regardless of the type of the variation of foF2, and that the features of the variation 
are the increase at night and small decrease in daytime. Since the height variation 
is independent of the type of the fof2 variation, the classification of the F2 distur- 
bances is done only for the f.'2 variations. 

The variations with season, sunspot or geomagnetic activities of the F2 distur- 
bances are shown in Fig. 4 for the four stations. In this figure the daily mean values 
(AfpF'2) of 4f,F2 on the first and second storm days are plotted against the sum of three 
hourly Kp indices of the day and against the sunspot number in the month when the 
storm occurs. The negative value of ‘AfpF2 corresponds to the negative disturbance 
and the positive value to the positive disturbance. The circles represent the summer, 
filled circles the winter and the crosses the equinox. From the figure the following 
points are recognized. 

(1) In-high latitudes the negative disturbance seems to occur predominantly in 


three seasons. 
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Fig. 4, (c) Fig. 4, (d) 
Eig. 4 (a)-(d) The Distribution of 4f,F2 (daily mean of 4f,F2) against SKp and against the 
sunspot number in the month when the disturbance occurs. 

(2) In middle latitudes the negative disutrbance seems to occur predominantly 
in summer or equinox and positive disturbance in winter. For the small sumspot or 
geomagnetic activities the positive disturbance tends to occur. 

(3) In low latitudés or in the equatorial zone the positive disturbance occurs 
predominantly in three seasons. 

(4) The amplitudes of the deviations in the two types increase with the incresing 
sunspot or geomagnetic activities. 

In other words, the F2 disturbances in middle latitudes have two types, one of 
which is a representative disturbance in high latitudes and the other is one in low 
latitudes. In the following sections we show, using the data in middle latitudes, that 
the individual disturbances of the two types can be explained by the drift theory. 


3. Method of Treatment of Geomagnetic Data and Theoretical Considerations 


We assume that the electron drift is considered to be caused by an electric field 
deduced from the disturbance-daily variation of the geomagnetic field and that this 
electric field exists throughout in the F and E regions. Since in the present case the 
individual state of the F2 disturbance is to be studied, the corresponding disturbence- 
daily variation of the geomagnetic field should be derived from one geomagnetic storm. 
However, it ‘is difficult to make such a deduction in one station, so that we adopt a 
following method for the sake of convenience. That is to say, the values of the hori- 
zontal intensity (H) and the declination (D) of the geomagnetic field for the local time 
nearest to that of the commencement of the geomagnetic storm, are taken and from 
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these values the mean values of.H and D on the five international quiet days corres- 
ponding to that hour are respectively subtracted. The values thus obtained are the 
deviative parts at hour zero of the storm time. The process is followed for succeeding 
60 hours (2 days and 12 hours). | They are denoted respectively as 4H and 4D. Next 
the overlapped mean of 24 hours of them are calculated for the former 48 hours, 
assuming that the deviative values before the commencement of the geomagnetic storm 
are zero. Then the variations can be said to correspond to the Ds¢ variation of the 
average state. When these variations are subtracted respectively from 4H and 4D, 
the parts corresponding to the Ds variations of the average state can be derived. They 
are the requested disturbance-daily variations of the geomagnetic field. 

When the disturbance-daily variation Ds is derived, the electric field and then the 
velocity of the electron drift can be derived, knowing the electrical conductivity of the 
ionosphere. When an anisotropic electrical conductivity is taken into account and it 
is assumed that the winds in the E and F regions have the same directions and 
magnitudes and that the geomagnetic field is produced by the dipole field, the vertical 
velocity of the electron in the F region is given by [20] 


fa ee eee gee sin 0 
o= lie +E brass cos? 6) (1) 


where J, and J, represent the total electric currents respectively in the x (south) and 
y (east) directions flowing in the E and F regions, 6 the geomagoetic co-latitude of the 
station, F the total intensity of the geomagnetic field and 


A A 
Kay= yo Kap, 


> 74 = eee ap RuRy, s ( 2 ) 


where k= ( O'xx AZ+ [ores dz, etc., (33) 


~E ry 


and ozz’ etc. are the elements of the electrical conductivity tensor o’ ae qh ‘Ene 
integration is carried out for the E and F regions. J, and J, are related approximately 
to the geomagnetic components in y and x directions, that is 4Y and 4X, in such a 
way that 
22 Go /fAy, 2x Iy= —fAX, (4) 

where f is constant and taken as unity. Jf 4X, 4Y are the values deduced from the 
Ds variations and K.,, Ky are the values on the distubed day, the v equal to the vg 
on the disturbed day. If the station situates on the magnetic equator the expression 
of the vertical velocity is simplified as follows, 
=A 5 Lema (5) 
where K is the electrical conductivity of the ionosphere on the magnetic equator. 

In the course of the study it is necessary to know the daily variation of the F2 


Vv 


region on the quiet day. The vertical drift of the electron on the quiet day is con- 


sidered to affect the quiet daily variation of the F2 region. The vertical velocity of 


the electron (v;) on the quiet day is obtained by the same expression as (1) or (5), 
replacing the electrical conductivity and the variations of the geomognetic field in their 
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expressions by the values on the quiet day. The values of 4X and 4Y are those 
derived from the mean value of the S, variations of the five international quiet days 
in the month when the storm occurs. 

If the effect of the horizontal velocity and the height-gradient of the vertical 
velocity are not taken into account, the variation-of the electron density is given by 


On 24 (t2)—B (2) 0) (6) 


where represents the electron density, v the vertical velocity of the electron, q and 8 
the production and effective attachment rates of the electron respectively, ¢ the time 
and z the height in the unit of the scale height. g is taken as qexp (1-z-e~* sec %), 
where q is constant and x is the solar zenith distance. The height distribution of Bp 
used for calculations is shown in Fig. 5. The equation (6) is solved by the numerical 
method firstly attempted by Millington [21]. Then the equation is replaced by the 
following forms i 


009g =F le |E ly a 
2am a 
sige v 


where t=1.37x10!¢, 1/oo= 1.37 x 1048, 1/%o 
=Bo/do, Y=N/M, B=Bo\€\z (Bo is constant) and 
|F|,=exp(1—z—e-* secy). If the velocity v 
equals to v, the solution of the equation 

| represents the daily variation on the quiet day 
and if v equals to vg+v_ the solution repre- 
sents the daily variation on the disturbed day. 
When the former variation is subtracted from 
the latter the residuum equals to the deviation 

~ from the normal variation. 
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Fable 1 p(t) on the quiet day 
local time | | 2 | 4 | 6 | 8 | 10 | 12 | 14 | 16 | 18 | 20 | 2 
coud “| 0.48 | 0.42 | 0.42 1.2 | 2.4 4.0 | 4.6). 4.0. 2.4 oe 1.0 | 0.6 
ey a P| eee be Bat toe BS ee |! | | | 
equinox | 0.43 | 0.41 | 041 wA.d-4 (2.2) 3.6 | 4.2 | 3.6 ee? “1.2 | 0.8 ‘ 0.5 
winter | 0.40 0.35 0.35 1.0 2.0 See 3.8 pe | 2200 1.0] 0.6 0.5 


Table 2 p(t) on the disturbed day 


local time o | 2 ) 4 | 6 | 8 | 10 | 12 | 14 | 16 | 18 | 20 | 22 
summer | eps visa) 32] 3.6| 4.0] 4.6] 4.0| 3.6] 3.2| 3.2] 3.0 
| equinox Par) 2.7 | 3.0.9.0 [be.2i)n9.6 |» 4:2) 8.618 3.34310 |/n3.0 oe 
| winter [oe irr sea Fpet7 2.9 | 'Bt3"|"3:8 3.3/2.9 2.9) 2.7| 2.8 


for three stations for the sake of simplicity. The daily variation factor p(t) varies with 
season. It is shown in Table 1 for the quiet day and in Table 2 for the disturbed day. 

In Table 2 we assume that the electrical cunductivity at night is comparable 
order of magnitude with that in daytime, though it is doubtful whether the assumption 
is correct or not. This assumption is derived from the fact that the geomagnetic 
variation on the disturbed day at night is equivalent to that in daytime. 

The method of the calculation of the deviative variations is stated in the section 
3. The results for some ionospheric disturbaces are shown in Fig. 6 (a)-(n)., In these 
figures the part A represents the observed variations of 4foF2 and Ah’F2 (or 4h;F 2), 
B, the corresponding calculated values, C, the disturbance-daily variations of 4H and 


Observed 


24 Bane () 24 hours \e boy 
Stormlime Storm Time ’ 
(a) Washington, 0’, Sep. 7, 1952 (b) Washington, 16%, Feb. 23, 1952 


0 24 hours 48 0 24 hours 
Storm Time Storm Time 
(c) Washington, 22”, Jan. 4, 1952 (d) Washington, 16”, Jun. 3, 1949 


éL 
4 
2 
° Mskecs ry 
StormTime Storm Time 
(e) A hap ay olin 1”, May 12, 1949 (f{) Watheroo, 3%, Dec. 16, 1944 


‘Fig. 6, (a)-(n) Observed and calculated deviations of individual negative disturbance. 
(A) represents the observed variations of 4f,F2 (dotted line, right ordinate), 4h’ F2 (or AhpF2 for 
kokubunji, full line, left ordinate) and the corrected variation of Ah’F2 (full line with circle) 
for the F2 disturbance which begins at the time shown below. The time is the nearest local 
time before the beginning of the geomagnetic storm and this time is taken as zero of the 
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(g) Watheroo, 10*, May 1, 1943 (h) Watherffo, 1”, Dec. 21, 1942 


48 
() 2 hours a 2 hours 
Storm Time StormTime 


(i) Watheroo, 20", Nov.23, 1942 (j) Watheroo, 15", Apr. 24, 1941 


storm time. (B) represents the calculated variations of the deviations of the maximum electron 
density (dotted line, right ordinate) and the height (full line), (C) the observed geomagnetic dis- 
turbance-daily variatons of 4H (full line, left ordinate) and JAD (dotted line) and (D) the foF2 and 


Storm Time 


(k) Watheroo, 17”, Mar. 28, 1941 


Saito eta eS 
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some examples are not obtained in our calculations. 

As shown in the section a it is evident that the large increase of the height in 
daytime is considered to be the retardation of the radio wave, so that the correction 
should be made for these values. The correction in Fig. 6 is made based on the theore- 
tical curve shown in Fig. 3 and the corrected variations are shown by a full line with 
circle. Then it is clear that the cal¢ulated variations of the height are consistent with 
the corrected variations. 

Conclusively it can be said that the F2 negative disturbances are well explained 
by the drift theory. 


5. F2 Positive Disturbances 


At first we must take care of the fact that two types of apparent positive dis- 
turbance may take place. One is a type where the fluctuation of the solar radiation 
responsible for formation of the £2 region leads to the apparent positive deviation. 
The other is a type seen when many days of the negative disturbance exist in a month, 
for the normal variation is mistaken as a positive disturbance, since the median value 
of the month becomes lower than the normal value. In the following the above appa- 
rent disturbances are excluded. 

Now, calculations of the deviations of the maximum electron density and the height 
for some positive disturbances are carried out by the same process as those for the nega- 
tive disturbance. The stations and periods of the data used are the same as those in 
the negative disturbance, and the electrical conductivities are those shown in (8). The 
results are shown in Fig. 7 (a)-(d). In these figures the parts A: to C correspond to 
those in Fig. 6. From figure 7 it is found that the calculated variation of the. height 
accords well with the observed, while the calculated variation of the electron density 
does not show sufficient coincidence, though the accordance seems to be good for the 
first day. 

It is important to know that the increase of the electron density is caused by 

the vertical drift of the electron. This is, as shown in the section 7, due to the mecha- 
nism that the drift velocity on the quiet day is reduced by the drift velocity additively 
generated on the disturbed day, so that the region tends to cease to move (the electron 
density on the static region is larger than that of the moving region). The insufficient 
coincidence in Fig. 7 is perhaps due to the fact that the daily variation of the drift 
velocity deduced from the S, variation on the quiet day is somewhat different from the 
actual state in the F2 region and the machanism of the increment does not become 
effective. The phase and the magnitude of the drift velocity depend on the values of. 
K,, and K,, and are modified by the induced field by the wind in the F2 region. 
Especially, the effect of the wind in the F2 region may give rise to the different phase 
or magnitude of the electric field in that region from those in the & region. (This 


- possibility is larger in winter than in summer, since the electron density in the former 


season is larger than that in the latter [29], [30].) We must also pay an attention to 
the fact that the wind from the observation in winter varies diurnally, while the 
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(a) Watheroo, 1”, May 18, 1943 


Storm Time 


(c) Kokubunji, 9%, Jul. 20, 1952 


is taken as zero of the storm time. 


observed geomagnetic disturbance-daily variations of JH (full line, left o 


line, right ordinate). 


(b) Watheroo, 2", May 24, 1940 


Storm Time 


(d) Kokubunji, 1%, Feb. 22, 1951 

Fig. 7, (a)-(d) Observed and calculated deviations of individual positive disturbance. 
(A) represents the observed variation of 4f,F 2 (dotted line, right ordinate) and 4h’ F2 (or AhpF2 
for kokubunji, full line, left ordinate) for the F2 disturbance which begins at the time shown below. 
The time is the nearest local time before the beginning of the geomagnetic storm and this time 


(B) represents the calculated variations of the deviations of 
the maximum electron density (dotted line, right ordinate) and the height (full line) and (C) the 


rdinate) and 4D (dotted 


S, variation (and also v,) in that season behaves predominantly semi-diurnally [27]. 


Disturbances in the Ionospheric F2 Region Associated 15 


6. F2 Disturbances in Equatorial Zone 


For the comparison with the results in the sections 4 and 5 and for the later 
discussions, we show the observed and calculated results of the 2 disturbances in the 
equatorial zone. These results aye published in another paper [9], hence the brief des- 
cription is given. The study is made for the data at Huancayo during the years 1938 
-1939. The electrical conductivity is taken as 


K=1.73X107* p@) emu (9) 


both for the quiet and disturbed days. p(¢) in (9) has no seasonal variation and is shown 
in Table 3. 


Table 3 p(t) in the equatorial zone 


local time | # | 2 | 4 | 6 8 10 | 12 | 14 | 16 | 18 | 20 | 2 


p(t) | 0.27 | 0.21 | 0.21 | 0.58 | 1.36 | 2.20 | 2.57 | 2.20 | 1.36 | 0.58 | 0.21 | 0°21 


Some examples of the calculated and observed deviations for individual state are shown 
in Fig. 8 (a)-(c). In this figure parts A and B correspond to those in Fig. 6 and 7, and C 
represents the disturbance-daily variation of 4H (~4X). 


Most of F2 disturbances in the equatorial zone are the positive disturbances and 
it is clear that these disturbances are well accounted for in term of the vertical electron 
drift both for the electron density and the height and that above results are due to a 
fact that the velocities of the drift on the quiet and disturbed days are set up with a 
phase difference by about 180°. 


7. Variation of F2 Disturbances with Season and Latitudes 


In the section 2 the observed variations of the F2 disturbances with season and 


Observed 


Storm Time Storm Time 
(a) Huancayo, 20", Jul. 2; 1939 (b) Huancayo, 19”, Jan. 4, 1939 


16 


(c) Huancayo, 16”, Feb. 13, 1938 


Fig. 8, (a)-(c) Observed and calculated devia- 
tions for the F2 disturbance in the equatorial 
zone. (A) represents the observed variatons 
of Af, 2 (dotted line, right ordinate) and Jh’ 


F2 (full line, left ordinate) for the disturbance 
which begins at the time shown below. The 
time is the nearest local time before the 
beginning of the geomagnetic storm and this 
time is taken as zero of the storm time. (B) 


repressents the calculated variations of cevia- 
tions of maximum electron density (dotted 
line, right ordinate) and the height (full line) 
and (C) the observed gomagnetic disturbance- 


ce} 24 hours 
Storm Time 


48 daily variatin of 4H. 


latitude are revealed and in the sections 4,5 and 6 they are accounted for as an effect of 


the vertical electron drift for individual cases. 


In this section we discuss physically why 


the two types of disturbance occur, why the occurrence of them depends on the season 


and latitude and furthermore 
why the amplitude of the devia- 
tive value depends on the 
magnetic or sunspot activities. 
Fig 9 shows the mean 
variation of the velocity v, on 
the-quiet day (dotted line) and 
the mean variation of the veloci- 
ty vg on the disturbed day (full 
line) in local summer (left hand 
side) and winter, at Washington, 
Watheroo, Kokubunji, 
yearly mean vaviation of those 


and 


velocities at Huancayo. From 
this figure it is easy to find that 
the velocity vg in each station 
has the same phase (with the 
maximum downward velocity 
near noon). The velocity vg in 
higher latitudes than Washing- 
ton will have the same phase, 
inferring from the analysis of 


Sp by Vestine et al [22]. On 


Summer Washington, Winter 


ed 


0 2 4 6 8B 10 [2 [4 16 18 20 22 24 
Local Time 


Fig. 9 Mean daily variations of the vertical drift velocities 

vq (dotted line) and vg (full line) for #2 disturbances used in 

the calculations. The variations in summer are shown in the 

left-hand side and those in winter and yearly mean at 
Huaneayo in the right-hand side. 
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the other hand, the velocities vw, in the four stations differ depending on the season 
and latitude. That is to say, in summer the velocity has the phase whose maximum 
downward value comes near noof at Washington, changes the phase with latitude and 
has a phase whose maximum upward value comes near noon at Huancayo, while in 
winter in middle latitude v, varies from 0" to 14” analogously to that at Huancayo. 

Now the maximum electron density in the moving region is lesser than that in 
the static region (no movement). This is well known as the geomagnetic distortion in 
the magnetic equator or in middle latitudes [23] [24] [25], Therefore if the velocity va 
is set up in such a way that vg+v, tends to zero, that is to say, that the v, differs 
from vg by 180° when their amplitudes are about the same, the moving region tends 
to approach to the static region and the maximum electron density increses. This 
representative example is the positive disturbance at Huancayo. Analogous circums- 
tances take place partially in winter in middle latitudes as shown in Fig. 9, so that 
the positive disturbance also takes place in those latitudes in winter. 

On the contrary, if the daily variation of vz has the same phase as that of v, the 
reduction of the maximum electron density is promoted. In this case the negative dis- 
turbance takes place. The representative example is that at Washington in summer 
or in equinox. For the negative disturbances in Watheroo and Kokubunji in summer, 
above principle may be applied, though the coincidence of the phases of both velocities 
is not sufficient. Since the verocity vg is deduced from the geomagneitc S, variation, 
which originate from the dynamo current in the ionophere, the latitudinal and seasonal 
variations of the S, daily variation leads to the variation of v, and to the daily variation 
of the F2 region. According to the discussions above mentioned, in the higher latitudes 
than the center of the S, current system the negtive disturbance may occur, since the 
daily variation of the deduced vg may coincide with that of vg, and in lower latitudes 
than the center the positive disturbance is easy to occur. In summer the center of the 
S, current system situates in fairly low latitudes, about in 30°, therefore in higher 
station than 30° the negative disturbance may occur, while in winter the center seems 
to situate near 40°, therefore below this latitude the positive disutrbance may predomi- 
nantly occur. (though actually we must take into consideration the magnitudes of the 


velocities v; and vq as described below). 
Consequently, the dependency on the latitude and season of the occurrence of 


the negative or the positive disturbances is ascribed to the variations of the drift 
velocity vz on the quiet day with latitude and season, in other words, to the effect of 
the S, dynamo current in the ionosphere. 

The dependency on the geomagnatic or sunspot activities of the amplitude of the 
deviation and of two kinds of disturbances is explained in the following way. With 
the increasing geomagnetic activity the disturbance-daily variation, from which the 
velocity vg is derived, increases. In this case the magnitude of vg becomes to be 
negligible compared with that of vg, so that the resultant velocity of vg and vg romoves 
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the seasonal variation and sometimes the latitudinal variation. After all, the drift 
velocity in the F2 region in the strom time becomes apparently vz alone and the decrease 
of the maximum electron density results in every latitudes and seasons. The larger the 
geomagnetic activity, the larger is this tendency. Thus the negative disturbance is easy 
to occur in large geomagnetie activity. 

There seems to be two reasons for the fact that in larger sunspot activity the 
larger (negative) deviation of the fpF'2 is caused. One of them is that in larger sunspot 
activity the fF 2 itself becomes larger, therefore the decrease of the fF2 must be 
large compared with that in low value of the fo'2, even if the decrease of the electron 
density has the same percentage. The other is that in larger sunspot activity the larger 
geomagnetic variation is easy to accur, so that the velocity vz becomes large compared 


with v, and thus the negative bisturbance tends to occur as mentioned above. 


8. F2 Disturbances in High Latitudes 


The theoretical considerations of the F2 disturbances at the stations in higher 
latitudes than Washington are not given, because of lack of the hourly values of the 
ionospheric or/and geomgnetic data in those stations. However, Meek [12] studied the 
ionospheric disturbances statistically in the high latitudes includuig the auroral zone 
(Canada and Alaska) for their individual states. He pointed out that the most remarkable 
feature of the F2 disturbances in these latitude is the depression of the f>F2, whose 
maximum comes at noon (somtimes the second maximum depression at about 18"). Thus 
the F2 disturbances in high latitudes are predominantly the negative disturbances. This 
is also known from Fig. 1 and 4 in the section 1. Meek also pointed out that the negative 
variation continues for several days and that the greatest depression of the f,.F'2 takes 
place in the southern part of the auroral zone, 60° to 65° geomagnetic latitudes, and 
not in the auroral zone itself. Since he did not discuss there the relation between the 
individual F2 disturbance and the corresponding geomagnetic variation, we cannot 


know the individual geomagnetic variation. But the general characteristics of the Sp 
variation in high latitudes (below the auroral zone) are similar to those in middle latitudes, 


as shown by Vestine et al. Hence the main feature of the F2 disturbance in high 


latitudes are considered to be accounted for as an effect of the vertical electron drift. 


9. Discussions 


(i) The effect of the electron production plays an important role on the deter- 
mination of the maximum electron density and the height both on the disturbed and 
quiet days. However, since the production has a complicated distribution with height 
and varies with time, it is sometimes neglected or elimibagted in the process of the 
calculation of the above values for the sake of simplicity. But this simplified method 
does not show the actual amounts, especially, for the height variation. In the sections 
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2, and 4 it is shown that the increment of the height in daytime on the disturbed 
day is not obtained. This is understood by knowing that the effect of the electron pro- 
duction is predominant. That is to say, the direction of the resultant velocity on the 
disturbed day is upwards in the early morning or sometimes throughout the morning, 
so that the electron mass of the Maximum density moves upwards approximately with 
the drift velocity when the electron production is absent. But when the production 
begins, the accumulation of the electron at lower height than the height of the initial 
maximum density, increases abruptly even the electron drift is going on, and the 
density there becomes larger than the initial maximum. Then the density and the 
height to be measured are those for the lower part. This displacement of the obserbed 
height occurs in a few hours after the sunrise. This circumstance will be the same 
for the drift of the electron by any causes other than the electric field. Consequently, 
the value of the deviations of the maximum density and the height can be obtained 
only by solving the continuity equation numerically or analytically which includes the 


electron production. 
(ii) It is found that the increase in height of the negative disturbance in 


daytime is not the actual variation and that as far as the height is concerned the 
daily variations of the deviation in most F2 disturbances have the same characteristic 
regardless of the latitude and season. Hence the statistical results for the height varia- 
tion which have been previously obtained by some authors, should be reexamined because 
of the correction of the height variation is not made. It has been considered up to date 
that since the increase of the height in many disturbances is accompanied by the 
decrease of the electron density and vice versa, the inverse relation between them may 
be the basic fact in the F2 disturbance. However, this concept no longer holds. 

(iii) In the statistical study of the F2 disturbance the deviation is often divided 
into two parts, Ds and Dst, by analogy to the analysis of the deviation of the geo- 
magneyic field. However, the separation in the ionospheric case seems to be brought 
as results of the technical treatment of the data and the physical meaning is some- 
what different from that of the geomagnetic case. According to our calculations, the 
vertical electron drift which has the diurnal variation causes the variations of the 
electron density and the height which contribute to the Ds and Dst variations. On 
the other hand, the electron drift which has a long period causes a variation which 
contributes to the two parts. These facts show that the variation in the ionospheric 
case differs from the electro-magnetic induction as seen in the geomagnetic field. The 
reason why the variation in the ionosphere is complicated is due to the fact that the 
variation depends on some factors in the continuity equation. From the above con- 
- siderations it is premature to consider that the Ds and Dst parts of the ionospheric 
-yariation are due to the causes which have respectively the characteristics like Ds and 


Dst. 
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(iv) Besides the phase difference of the drift velocities, the mechanism of the 


increment of the electron density are studied. They. are :— 


(a) The height gradient of the drift velocity 
(b) The temperature effect 
(c) Theelectron production by impinging ‘particles 


Of these considerations, (c) may be neglected for the disturbances in middle latitudes. 
In winter the height of the F2 region on the quiet day lies in the lower level than in 
other seasons, so that the effect (b) may take pace, because the temperature below 
300 km altitude has the gradient with the height [26]. This effect is calculated numeri- 
cally with an assumption that the temperature is 600 K at 200km and the linear increas- 
ing rate of the temperature is 5°K/km. It is found that the effect of the temperature 
gradient is negligible. For the effect (a) the calculations are made for the case of the 
increasing velocity upwards. ‘According to Briggs and Spencer [27], the horizontal 
velocity of the wind in the F2 region has a height gradient, of about 1m/s /km. As 
the wind in the F2 region is considered to be the motion of the charged particles due 
to the electric and geomagnetic fields, the motion in the vertical direction may has also 
a height gradient. So that the velocity seems to increase with the increasing height. 


The effects are calculated for the following three cases. 


Ql) vn 21? 


| (2) v0 = 
~ » 3) 00 G1) pee 


where uv» is the velocity at the reference level (g=0 for 2 =2-for (2), and ae (3)) 
and in each case the ieryelociy below the datum level is taken to be equal to the e velocity at 
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turbance can be explained by the drift theory. From the above results and the 
successful explanation of the “‘F 2 disturbance in the equatorial zone, it can be said 
that the negative disturbances Js caused most remarkably when the phase of the daily 
variations of the velocities v; and v, are about the same and the positive disturbance 
caused when they are different by about 180°, and that the seasonal and latitudinal 
variations of the F2 disturbance are attributed to the seasonal and latitudinal variations 


of the vz, in other words, to the S; dynamo current from which v, is deduced. 
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Abstract 


Palaeomagnetic studies have been conducted of the volcanic rocks in the 
North-Izu and Hakone volcanic region, Japan, where complete succession of lavas 
has been determined by one of us (H.K.). By sampling 4~7 oriented rock speci- 
mens at each of 57 sites, the period from the very beginning of the Pleistocene 
to Holocene has been covered. The maximum time interval between two con- 
secutive samples may probably be not more than several tens of thousand years 
except that between two samples of middle to younger Pleistocene when the 
volcanic activity did not occur within the region concerned. Care was taken not 
to use the rock samples of which natural remanent magnetization may have 
suffered from any significant disturbances, geolgically, chemically, magnetically 
or otherwise. Selection of proper-samples was performed according to the criteria 
for the stability of remanent magnetization proposed by us previously (Journ. 
Geomag. Geoelec., 6, No. 4, 1955). The major findings obtained through the 
present study are: 1) During the whole Quaternary age, the geomagnetic centred 
dipole has been fluctuating around an axis of which north pole wandered from 
72°N, 86°E to 81°N, 32°W. 2) The direction of polarization of the centred dipole 
was reversed at a time in the earliest Quaternary, namely during the middle period 


of the formation of the Usami volcano. 


1. Introduction 


Recent investigations on palaeomagnetism have brought about many findings 
of great geophysical importance, especially with regards to the long controversial 
hypotheses of the polar wandering, the continental drift and the reversals of the geo- 
‘magnetic field in the past geological time. Excellent summaries on these works are 
found in the publications of S.K. Runcorn [1], P.MS. Blackett [2] and E. Irving [3]. 
Many authors, on the basis of their own measurements and others’, have suggested 

that the earth’s magnetic field had undergone repeated total reversals during the geo- 
‘ logical time. Such suggestions, though plausible, would be completely justified as 
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physical truth only when the contemporeinity of reversely magnetized rocks from 
various localities of the world could be proved with sufficient certainty for each sug- 
gested reversals. Such absolute datings of rocks are, however, almost inaccessible at 
the present time. In general, the ages of rocks are determined by geological means 
more precisely for the younger strata. Therefore, since 1954 we have engaged our- 
selves with a palaeomagnetic research on a Quaternary volcanic rocks of North-Izu and 


Hakone volcanic region (Fig. 1), where the complete succession of lavas had been deter- 


\ 
North Izu 
and Hakone Regi 


Fig. 1 Map showing the locality of the North-Izu and 
ab Hakone volcanic region. 
mined by one of us (H.K.) [4], with the hope that our data would possibly be compared 
with the data for other parts of the world such as Iceland and France, where reversely 
magnetized rocks of early Quaternary were reported by J. Hospers (5h ay -Einarsson 
and Th. Sigurgeirsson [6] and A. Roche [7]. It must here be noted. that Y. Kato in 
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Table 1. Stratigraphic succession of rocks of North-Izu and Hakone region and 
the estimated period of their eruption. 
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was covered by an open sea where sediments containing mulluscan fossils were deposit- 
ed (Zy6d formation). 

The oldest of the Quaternary volcanoes is Usami which was erupted from a 
centre near the present village of Usami, north of Ito. It is a stratovolcano of pyro- 
xene andesite (designated as UV in the tables). Some of the lavas which flowed south- 
westward covered the Zy6 formation. 

Slightly later the centre of eruption moved about 10 km northward where another 
stratovolcano Taga was born. From this centre lavas and pyroclastic materials of 
pyroxene andesite were erupted in the earlier stage of activity, whereas fluidal lavas of 
olivine-pyroxene basalt, almost to the exclusion of pyroclastic materials, were extruded 
in the later stage (all designated as TV). 

After a short period of erosion lava domes of dacite were erupted near the centre 


of Taga and also on its northwestern flank. f 


Then the centre of activity migrated again about 10km north-northwestward, 
namely to a point near the present town of Yugawara, from which lavas ‘and pyro- 
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clastic materials of pyroxene andesite were erupted to form Yugawara Volcano (YV). 

Slightly later the activity of’ Hakone Volcano started from a centre situated about 
10km north-northwest of the town of Yugawara. The activity of this volcano took 
place in three distinct stages, the two contiguous stages having been separated from 
each other by periods of caldera formation and erosion. 

During the first stage of activity a stratovolcano about 2700m high was built. 
This eruption was caused at the outset by a magma of olivine-pyroxene basalt and 
later by that of pyroxene andesite. These rocks are called the Old Somma lavas (OS). 
Then the top of the cone subsided to form caldera measuring 11 km in NS diameter and 
7km in EW diameter. This was followed by a long period of erosion during which 
more than 800m of the Old Somma lavas and the underlying Tertiary rocks in the 


~ eastern caldera wall were removed. 


The second stage of activity was caused by magmas of salic pyroxene andesite 
and pyroxene dacite (the Ycung Somma lavas, YS), resulting in a gentle-slcred shield 
volcano within the caldera already formed. At the close of this activity a tremendcus 
amount of pyroxene dacite magma was drained from the reservoir and produced pumice 
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Fig. 4, 1 Fig. 4, j 
Fig. 4,a Sample No. 26, andesite of Umenokidaira volcano. 
Fig. 4,b Sample No. 9, Ov, andesite of Old Somma of Hakone volcano. 
Fig. 4,c Sample No. 57, Oj, andesite of Old Somma of Hakone volcano. 
Fig. 4,d Sample No. 16, D,, pyroxene dacite at Akanezaki. 
Fig. 4,e Sample No. 17, TV;, basalt of Taga volcano. ; 
Fig. 4,f Sample No. 15, TV;, andesite tuff breccia of Taga volcano. 
Fig. 4, g Sample Nos.18 and 45, UV, reversely magnetized andesite of Usami volcano. 
Fig. 4,h Sample Nos. 20, 21 and 22, UV, normally magnetized andesite of Usami volcano, 
Fig. 4,i Sample No. 3, CC;, andesite of a central cone of Hakone volcano: example of 


poorer uniformity in direction of natural remanent magnetism. 
. _ Fig. 4,j Sample Nos. 44 and 46, fragmentary part of andesite of Usami volcano : direction 
of natural remanent magnetism is widely scattered. : 
Fig. 4 Examples of direction of natural remanent magnetism in Schmidt projection. 
“Full circle... lower hemisphere 
Hollow circle .. upper hemisphere 
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flows. This eruption was succeeded by collapse of the centre of the shield volcano, 
resulting in a caldera approximating the earlier one in size and location. 

After a long period of erosion which gave rise to gorges 300m deep cutting 
through the eastern flank of the shield, the activity of the third stage started. Viscous 
magma of pyroxene andesite was extruded from vents arranged on or near an NW-SE 
line passing through the centre of the caldera. Six lava domes and one stratovolcano, 
which are called the Central Cones (CC), were formed by this activity. 

The materials erupted at various stages during the growth of the Hokone Volcano 
are interbedded with fossiliferous sediments exposed in a region to its northeast. From 
stratigraphical and palaeontological data obtained from this region, it is infered that 
the birth of the volcano was probably middle Pleistocene. The materials of the Central 
Cones is interfingered with younger Pleistocene lapilli and ash from Huzi (Fuji) Volcano 
lying to the northwest of Hakone. It follows that the activity of Hakone lasted for 
the later half of the Pleistocene. 

To the south of the Usami Volcano, eruption of fluidal basaltic lavas and basaltic 
scoriae took place during a period from younger Pleistocene to Holocene. These 
materials built up small lava platforms, shield volcanoes, stratovolcanoes, and pyroclastic 
cones within a area about 12 km in maximum diameter. These volcanoes are collective- 
ly called the Omuro-yama Volcano Group. 

In Table 1 the approximate thickness of the materials composing each unit of 
Pleistocene is given, together with the-period in 10? years during which each unit was 
formed. Each erosion period is roughly estimated in the following ways: 

The lavas and pyroclastic materials of Usami, Taga, Yugawara and the Old 
Somma of Hakone are similar to one another in character. It is assumed therefore 
that each of these units was accumulated during a period proportional to the thickness 
of its materials. Thus the earlier half of the Pleistocene (approximately 500,000 years) 
is divided into three periods according to the thickness of the rocks of Usami, Taga, 
and Yugawara, separated by short erosion intervals. The period of the Old Somma lava 
eruption is also estimated according to the thickness of the rocks relative to those of 
the older units. As the rocks of the Young Somma and Central Cones have quite 
different lithologic characters from those of the older units, the periods of their eruption 
may not be proportional to their thicknesses. Each period is estimated as 50,000 years 
regardless of their thickness. The erosion period between the eruptions of the Old Somma 
and Young Somma lavas and that between the eruptions of the latter and the Central 


Cones are roughly estimated as 200,000 and 50,000 years respectively on the basis of 
the volumes of materials removed. 


. 


3. Results of Measurements on Natural Remanent Magnetism 


Although each unit listed in Table | is composed of many layers of lavas and 
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pyroclastic materials, some of them are not suitable for the palaeomagnetic study, 
because they are either made up to fragments that had already been cooled through 
the Curie point before they came to settle at the present positions (pyroclastic rocks or 
fragmental parts of some lavas)’ or subjected to disturbances after the cooling. All 
samples were therefore collected from solid interior of undisturbed lava flows. One 
sedimentary rock is also included. These samples were selected so as to represent 
different horizons within each unit. On each outcrop 4 to 7 oriented samples were 
collected at points located more than 1m apart from one another. The localities of 
these outcrops are shown in Fig. 3. 

The direction and intensity of the natural remanent magnetism were measured 
for most samples by an astatic magnetometer, whereas those of a few samples with 
weak remanent magnetism were measured by a spinner type magnetometer. 

Some of the results of measurements are illustrated in Fig. 4, a, b,c, d,e, fools 
and j by the Schmidt projection. They illustrate the degree of concentration of mea- 
sured directon of the remanent magnetism of individual specimens at a site. Fig. 4, 
a, b,c, d,e, f, g and h are good examples, Fig. 4, i is the example of poorer concentration 
and Fig 4, j shows how the directons of remanent magnetism are scattered for a group 
of samples collected from a fragmental part of a lava flow. 

Table 2 lists the general 
results of measurements. The 
names of the geologic units and 
rocks refer to those in H. 
Kuno’s paper [4]. 

The values of direction 
and intensity of natural re- 
manent magnetism for each 
site and the error angles for 5% 
of direction were calculated 
from the measurements of 4 to 
7 samples by Fisher’s method 
[9]. The layers with cross 
marks are omitted ones on 
account of geologcal evidence 
as mentioned already. The 
direction of the geomagnetic 
centred dipole, listed in the 7th 
column of the Table 2, were 


Fig. 5 Directions of geomagnetic centred dipoles deduced 


calculated from the measured from the mean directions of natural remanent magnetism 

direction of the remanent at each outcrops: full circle....position of north pole 
hollow circle........position of south pole 

magnetism. They are plotted CONELE Marien sein ees present geographic pole 


in Fig. 5 in which the full and . 
hollow circles stand for the north and the south poles respectively. It is noted in 
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Table 2 and Fig. 5 that the direction of the geomagnetic field has always been the 
same with that of today during the whole Quaternary age, except for a short period 
when some of the lavas of the Usami volcano were erupted. During this period of 
early Quaternary the direction of geomagnetic field was nearly reversed to that of 


today. 


4. The Reliability Tests 


As we have discussed previously [10], the palaeomagnetic study by means of 
the natural remanent magnetism of an igneous rock is physically justified when and 
only when the direction of the natural remanent magnetism is ascertained to have 
been kept unchanged throughout the time after the rock had been formed. From the 
present knowledge of rock-magnetism, the above condition may be regarded as satisfied 
when the natural remanent magnetism in question is proved to be the stable 
thermo-remanent magnetism of a stable titanomagnetite phase, produced when the rock 
cooled through the Curie point at the time of formation. Quite recently, the existence 
of stable chemical remanent magnetization and piezoremanent magnetization has been 
reported [11] [12]. But, their stability has not thoroughly studied yet and, besides, they 
are mainly concerned with the remanent magnetization of sedimentary and metamorphic 
rocks which we are not dealing with here. We have proposed in the previous paper 
that a set of practical tests should be applied on the rock specimen, when one wants 
to use its natural remanent magnetism for palaeomagnetic purpose. The details of 
the tests being described in that paper [10], we shall only briefly explain them in items 
here. Since the following criteria are qualitative in nature, some quantitative con- 
ditions were tentatively adopted in the present study, as shown in items 1) and 4). 
Symbols ©, A and & mean the “ adoptable,” “somewhat dubious” and “ rejectable ” 


respectively. 
Notations to be used are as follows: 

HERG IES, 2c natural remanent magnetism at a temperature 7. 

TaD 568. thermo-remanent magnetism at a temperature 7. 

Le As oc wee room temperature. 

ACS, Se Curie point temperature. 

dS eg mode of thermal demagnetization of J, by heating up to T and 
cooling down to 7, in non-magnetic space. 

PTs es mode of thermal demagnetization of /r, by the same process as 
Needs 

ati a ee mode of demagnetization of J, by alternating magnetic field of 
which maximum intensity is H. 

J baie’ Portes mode of demagnetization of /7, by the same process as in (J,,)H. 

Uh Ta) Pah eae. saturation magnetization at a temperature 7. 


The reliability tests 
1) The ratio J,(T.)//re(T.) should be nearly equal to unity. 
Ou. . MOSS ii eho . 
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Bian ve O2<//re<05 Or'US <J,//7.<2.5 
Seer otherwise 
2) The modes of /,-T and J7,-£ should be similar. 
3) The modes of ( T,.)H and (J 7H should be similar. . 
4) Contained ferromagnetic mifieral should exclusively be a single phased titanium-poor 
titanomagnetite. 4 


Re ost ess favor T.>400°C 
Buian¥ ts = 400°C >T,>300°C , 
pT ee 300°C>T, 


5) J,(T) curve measured in vacuum should be reversible with temperature. 

The results of these tests are listed in Table 2 by the symbols. As these symbols 
show, there is no particular indication that the naturally reversely magnetized rocks 
have any properties which may cause the self-reversal phenomenon. Hence, we find 
no reason to doubt the reversal of the geomagnetic field at the time when these rocks 


uten 30 etext (a, arts § 


zy J ese 
se 


Mise tenrioe =k. 


e i 
SER Cs dP 
ae ola 
ae pe gt 


34 T. Nacara ef al. 


Ca 200 — 4 "CoO mis 200 7 hat 
ee eo “~~ Hin Oe. 


< 4 
4 
: ? 
iy xT 
‘ asm 3 as 
1-2 
x 4 4 3 
. 
a 2 4 7 
ma Pe na ima ? 7 
' - a Pen 
f f , 
o j z y. i 
. -. - £ ; 2 


' 
wy 
; 
n 
id 
rawr ae 
; 
, 


Palaeomagnetic Studies on a Quaternary Volcanic Region in Japan 5) 


36 ; T. Nacata ef al. 


Palaeomagnetic Studies on a Quaternary Volcanic Region in Japan ot 


Usr, . Jn 
, ' Ci 


aly A in Oe, 


(a) (5) 


Js(T) 
Jg(To) 


1.0 


(c ) Fig. 10 Reliability test of natural remanent 
magnetism: sample: No. 18, UV, andesite 
Test 1) Jn/JTc=2.4 
N 2) see Fig. 10, a 
Ne 3) see Fig. 10, b 
4) .see Fig. 10, c 
5) f= ol0-€ 
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were cooled through Curie point. Fig. 6, Fig. 7, Fig. 8, Fig. 9 and Fig. 10 are the 
examples of the tests performed. Fig. 10 shows the results of the tests for the reversely 


magnetized rock. 
5. Discussion and Conclusion 


The samples used in the present study cover the period throughout the Quater- 
nary age. Any pair of samples collected from nearest horizons represents a time inter- 
val which may not exceed several tens of thousand years except those representing the 
horizons immediately below and above the long erosion intervals in the history of the 
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Hakone Volcano. 
Examining the results obtained, it may reasonably be concluded from Fig. 5 that 
the sense of the direction of the geomagnetic field has remained the same as that of to- 
day throughout the Quaternary age, except for a short period during which a few 
successive lava flows of the Usami Volcano were extruded, and that the pole wandered 
about the neighbourhood of the present north with periods shorter than the duration of 
the volcanic activity which produced each unit. In order to show the more general ten- 
dency of the polar wandering during the Quaternary age, these rather rapid fluctuations 
are smoothed out and the value of the geomagnetic field representative of each 
period of volcanic activity is plotted in Fig. 11. It seems likely, judged from recent 
geomagnetic secular variation 
data and results of palaeomagne- So 
tic studies on historic and pre- : Geographycal Meridian 
historic baked earth that pseudo- 


ata the : Present Declination 
periodic fluctuations of the geo- gl 


magnetic pole of much shorter 2 e 
period and of less amplitude e 
are superposed on this smoothed ° 
curve of general tendency. The 
absolute time-scale taken for «© 
abscissa of Fig. 11 is obtained - 
as a mean value of the time 
range assigned to each unit given 
in Table 1. The value for the 
ordinate is the declination mea- 
sured from the present geogra- 
phical north and the inclination 
representative of each period 
of activity. They are obtained 
by averaging the values for the 
localities of the samples belong- 


ing to the individual periods’ 3 Mh er ee : 
& periods of Fig. 11 .,VariationSin the direction of geomagnetic field in 


volcanic activity. For the the North-Izu and Hakone region through the Quaternary 
Barcnk Golan tihe dunatiom of age, as deduced from palaeomagnetic evidence. 
activity was divided into three parts because of the occurrence of reversal of the 
geomagnetic field. Fig. 12 shows the movement of the geomagnetic north pole calcu- 
lated from the values in Fig. 11; for the reversal time the south pole is plotted by a 
hollow circle. Viewed from the north, the pole is thus estimated to have shifted 
-clockwise from 72°N, 86°E to 81°N, 32°W. The values shown in Fig. 11 and Fig. 12 
-are listed in Table 3. . non She ; ‘ 
In order to determine the duration through which the reversal lasted, we measured 
the direction of natural remanent magnetization for many samples of the Usami lavas 


Palaeomagnetic Studies,on q Quaternary 


> ees} Position of the present geomagnetic pole 


Fig. 12 Position of geomagnetic north (full circle) 
and south (hollow circle) pole in the Quaternary age 
as deduced from palaeomagnetic evidence: numbers 
attached to circles, as in Table 3, represent the geo- 
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(see Table 2). It is remarkable that 
reversal of natural remanent magne- 
tization is found exclusively in the 
lava flows occurring at the middle 
horizon of the lava sequence of the 
volcano. As the duration of the 
activity of this volcano is estimated 
as about 100,000 years, the length 
of the duration of reversal should be 
a fraction of it. 

It has not been possible, so 
far, to detect the possible change 
of direction of the transient geo- 
magnetic field, in spite of our nume- 
rous measurements for the layers of 
lava that had considered to form the 
boundary between normally and 
reversely magnetized ones. For 
example, along the road-cut imme- 


logic units. 
diately north of Usami Tunnel, 
Table 3. 

os Votoaie D, 1, error angle for ee pene fieal 

p=0.05 | p=0.10| p=0.50 coordinates.) 
8 OV -3° 44° 14° cial Bs 81°N 32°W 
7 CC -9 55 23 19 9 86 N 65 E 
6 YS -7 43 34 11 8 79 N 13 W 

Erosion 

5 OS -11 53 11 9 5 83 N 50 E 
; YV -3 46 - - - 82 N 30 W 
4 Vs -19 5S: 20 17 8 78 N 62 E 
8 UV; -33 59 15 13 6 70 N GIP E, 
2 UV, 146 -43 29 23 13 64S 145 W 
ib UV; -27 63 12 10 5 72 N 86 E 


northeast of the village of Usami, five successive lava flows are exposed without 
intervening layers of pyroclastic rocks, Complete absence of weathering products 


at the top of each flow strongly indicates 


that they were emplaced in rapid 


succession. All the flows but the upper-most one show the reverse direction of magne- 


tization. The degree of uniformity of remanent magnet 


ization of the upper-most layer 


is, as seen in Fig. 13, too poor to be adopted for the palaeomagnetic purpose, especially 


for such a critical use as a fossil of the geomagnetic field in a transient state from 


reverse to normal direction. The observed randomness in the direction of natural 
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remanent magnetization is regarded as due to the fragmentary nature of the layer. 

The present theoretical considera- 
1 tions, based on the self-exciting dynamo 
model of the geomagnetic field, on the 
other hand, suggest that the reversal of 
the geomagnetic field would have occurr- 
ed through vanishing and development of 
the intensity of the dipole field without 
accompanied by change in direction, 
rather than by the migration of the axis 
of the dipole relative to the axis of the 
earth’s rotation [13]. Our failure in 
tracing the intermediate directon of the 
geomagnetic field is, therefore, in agree- 
ment with these current views. We could 
also not find any systematic variation in 
Fig. 13 Direction of natural remanent magnetism the intensity"of remanent magnetization 
of the upper-most layer near Usami Tunnel: sam- Of normally and reversely magnetized 


ple No. 44, full circle....lower hemisphere 


; ded, from 
hollow circle....upper hemisphere. rocks. It may be concluded om the 


result obtained so far, that the change 
from the normality to reversal or vice versa took place fairly suddenly in comparison 
with the rate of accumulation of lava flows. More detailed investigations on the 
remanent magnetization of rocks emplaced at the reverse-normal boundary are now 
in progress. Evidence of sudden reversal of remanent magnetization has recently been 
reported by M. Almond, J. Clegg and J. C. Jaeger for Tasmanian Tuff of Triassic age 
[14]. 
Finally, it is of much interest to compare our data with those obtained in Iceland 
by J. Hospers [5] and T. Einarsson and Th. Sigurgeirsson [6], and in France by A. 
Roche [7]. Both data seem to be in good agreement with ours in that the last 
occurrence of the reversal of the geomagnetic field was in the early Quaternary age. 
But, as for the duration of the reversal, our data obtained so far seem to suggest that 
it was shorter than that suggested by the data in Iceland. 
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Abstract 


The geomagnetic secular variation during the period from 1950 to 1955 is 
estimated on the basis of annual mean values of the geomagnetic field at some 
seventy magnetic observatories distributed all over the earth. Judging from the 
spherical harmonic analysis of the variation, it seems likely that the magnetic 


moment of the earth’s dipole is still decreasing. 


1. Introduction 


The question of whether or not the magnetic moment of the earth’s dipole has 
recently begun to increase is of general geophysical interest and has been discussed 
by C. Gabier-Puertas [1], E.C. Bullard [2] and H. G. Macht [3]. On the basis of his 
extensive analysis [4], E.H. Vestine [5] showed that the rate of the decrease of the 
axial component of the earth’s magnetic moment becomes small after its maximum 
around 1920. On extrapolating his result, it has been suggested that the decreasing 
rate might become zero at some time around 1950 resulting in an increase of the 
moment after that epoch. 

Vestine’s analysis of the geomagnetic field at epochs of every ten years has also 
been the source of various studies related to the geomagnetic secular variation such 
as the westward drift, the electric current systems on the surface of the earth’s core, 
the shielding effect of the earth’s mantle and so on. In view of such geophysical 
significance of the geomagnetic secular variation, it will be of some use to make clear 
the recent status of the geomagnetic field, because some ten years have past since the 
last epoch of Vestine’s analysis. 

One of the writers (T.N.) has been appointed as chairman of the Special Com- 
mittee on Secular Variation and Palaeomagnetism, International Association of Geo- 
magnetism and Aeronomy since 1954. At his request, some fifty magnetic observatories 
kindly sent him the annual mean values of the three geomagnetic components during 
the period from 1945 to 1955. The writers can also add some additional data published 
elsewhere [6], [7], [8]. The total number of observatories whose data can be used to 
investigate the secular variation during the period from 1950 to 1955, amounts to es 
many as 68, so that the writers would here like to analyse the data in the customary 
way. The locations of these observatories are shown in Fig. 1. 


(42) 


Pope eT 


Geomagnetic Secular Variation. during the Period from 1950 to 1955 43 


o° 30° 60° 90° 20° + 150° 180 1580 = 120° e 


To 2 


22 


Fig. 1. The distribution of the magnetic observatories from which the data are supplied. 


Table 1. The changes in the three elements of the earth’s magnetic field at respective observatories. 
; (Unit : gammas) ; 
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(Table 1 continued) 


Observatory | Sener | Gorgiae | 4x av Az 
Manhay 50.°3 b7 124 215 107 
Prthonice 50.0 14.6 60 177 145 
Fiirstenfeldbruck 48.2 120m" 74 195 116 
Chambon-la-Forét 48 .0 2nd 147 201 59 
Nantes 47.3 358.4 174 203 72 
Yuzhno Sakhalinsk 46.9 142.7 103 —119 18 
Odessa 46.8 30.9 12 100 246 
Memambetsu 43.9 144.2 86 —42 —84 
Agincourt 43.8 280.7 ial at —144 
Vladivostok 43.1 131.9 83 —64 _ 
Dusheti 42.1 44.7 18 49 260 
Keles eel e 69.2 175 -19 277 
Istanbul-Kandilli 4nt | rs. gael 8 46 134 254 
Ebro 40.8 0.5 163 241 —2 
Coimbra 40.2 351.6 235 Zo —16 
Toledo 39.9 355.9 240 234 5, esd 
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Geomagnetic Secular Variation during the Period from 1950 to 1955 A5 
2. Data 


From the annual mean valties reported from each observatory, the changes in 
the north (X), east (Y) and vertical (Z) components during the period from 1950 to 
1955 are calculated as shown in Table I together with the geographical latitude and 
longitude of the observatories. Most of the values of X and Y are deduced from the 
data reported as those of horizontal intensity (H) and declination (D), while very few 
observatories reported the values of Z as combinations of H with inclination (J). For 
certain observatories, extrapolations or interpolations are necessary in order to calculate 
the values at the epochs required though the errors due to these procedures would 


seldom exceed a few gammas. 


3. Secular variation for the period 1950-1955 


The mean secular variation for the period 1950-1955 thus calculated is shown 
in three charts respectively for x. Y and Z as can be seen in Figs. 2a, 3a and 4a. 
In order to compare these charts to those of Vestine’s study for 1942.5, similar charts 
are reproduced from Vestine’s results as are also shown in Figs. 2b, 3b and 4b for 
respective components. 

On examining these figures, we see that the general mode of the secular varia- 
tion is not greatly different for both the epochs, though there can be seen some 
difference especially in Z. In Fig. 2a, we see three vortices for X centred at the North 
and South Atlantic Ocean and North India. The locations of these vortices are almost 
the same as those of Fig. 2b though the maximum intensities at their centres have 
changed within a range of a few tens of gammas per year. The negative centre 
which has been shown by Vestine in the southern part of the Indian Ocean seems to 
disappear. The distribution of Y is also nearly the same as that of the Vestine’s epoch 
as can be seen in Figs. 3a and 3b. The positive centre in North Africa and the negative 
ones near Madagascar and in South America still continue their activities, while the 
positive centre in the South Pacific seems likely to shift to the north somewhere around 
New Zealand. As for Zin Figs. 4a and Ab, the positive centre in the Middle East is still 
alive though the intensity becomes nearly one half. The negative in the Atlantic has 
moved to the north-west direction. We also see that the enormous positive extending 
from South Africa to the South Pacific around the Antarctica keeps its activity. Although 
the negative in the Indian Ocean seems to move to the south, no conclusive dis- 
cussion is possible because of the low accuracy of the data around this region as has 
been pointed out by Vestine. Cee: 

It is also interesting to compare Figs. 2a, 3a, and 4a to the non-dipole field for 
1945 which has been graphically shown by E.C. Bullard et al [9] by use of Vestine’s 
analysis. Most of the highs and lows of each component of the non-dipole field have 


been generally intensified during recent years. But the centres of the secular variation 


do not agree with these highs and lows. There is a marked tendency for the centres 
of the secular variation to occur at places to the west of the highs and lows of the 
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Fig. 2a. Secular vrriation for the period 1950-1955, north component, contour interval 10 gammas 
per year. 
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Fig. 2b. Secular variation for the period 1940-1945, north component, contour i 


nterval 10 gammas 
per year. (Reproduced from Vestine’s result.) , 
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Fig. 3a. Secular variation for the period 1950-1955, east component, contour interval 10 gammas 
per year. 
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or the period 1940-1945, east component, contour interval 10 gammas 
per year. (Reproduced from Vestine’ s result.) 


Fig. 3b. Secular variation f 
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Fig. 4a. Secular variation for the period 1950-1955, vertical component, contour interval 20 
gammas per year. 
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= 4b. Secular variation for the period 1940-1945, vertical component, contour interval 20 
gammas per year. (Reproduced from Vestine’ s result.) 
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non-dipole field. The fact can be interpreted as the westward drift of the geomagne- 
tic field as has been pointed out by a number of geomagneticians. However, no 
attempt to estimate the velocity of westward drift will be made in this paper, because 
the velocity is so small that such a short period as that covered by the present analysis 
would not give any significant results. 


4. Spherical harmonic analysis of the secular variation 


From Figs. 2a, 3a and 4a, the amounts of the secular variations for each 
component are read off at the intersection points of the meridians of every 30° and 
the parallel circles of every 15°. Using these values, spherical harmonic analyses are 
conducted in the customary way. The time-derivatives of the Gaussian coefficients 
can be calculated independently from the combination of X and Y and us the results 
being shown in Table II together with the previous analyses by many authors. 


Table Il. Spherical harmonic amplitudes of the secular variation for various epochs. 


P| ees, [| sees | RSS | _ventnetanee Sree 

| 1922-1885 | 1920-1902 | 1920-1902 | 1912.5 | 1922.5 | 1982.5 | 1942.5 verre 
93 | 20 42 0 25 28 23 9 9 6 
g:| -1 | -9 | 13 1 4 1 2 0 2 
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he Sir (S10 21s 10 ee ae 
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(Unit : gammas per year) 


As is seen in the table, the coefficients deduced from Ze agree fairly well with 
those from the combination of X and Y. Since the data concerned are supplied only 
from permanent magnetic observatories, the coefficients which are deduced from the 
vertical component are thought fairly reliable though they have been regarded as less 
accurate in the previous analyses. 
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5. Discussion 


Now we are in a position to compare the present result to that of the previous 
analyses. In an extensive investigation on all the spherical harmonic analyses available 
by 1945, however, P. Mauersberger [10] has emphasized the fact that no one analysis is 
strictly comparable with another. In fact, -the difference of the source of magnetic 
data, the method of computation and so on would result in some difficulties in com- 
parison of different analyses. As far as the relatively large coefficients are concerned, 
however, two analyses by different authors, for instance, the analysis by Vestine for 
1940-1945 and the present one, may be compared with one another in spite of the 
above difficulties. 

Since the mode of distribution is not greatly different as can be seen in Figs. 
2~4, the coefficients for both the analyses should agree fairly well with one another. 
In reality, we see in Table II that the relatively large coefficients such as 9°, 92, hi 
and h2 are in good agreement for both the analyses. Attention should be paid to the 
fact that g° is the same in both the analyses, so that the decreasing rate of the axial 
component of the earth’s dipole is to be the same for both the epoch. Even if possible 
errors of the coefficients are taken into account, it is of extreme doubt that the strength 
of the earth’s dipole began to increase some time after 1950. 


6. Conclusion 


The analysis of the secular variation of the earth’s magnetic field for 1950-1955 
shows that its general tendency is nearly the same as that for 1940-1945. It is most 
unlikely that the strength of the earth’s magnetic moment has been increasing some 
time after 1950. Instead, it seems still decreasing with a rate nearly equal to that for 
1940-1945. 

The writers are grateful to all the members of the magnetic observatories who 
kindly supplied magnetic data to the writers. The writers also wish to express their 
thanks to Miss E. Nakagawa who assisted them in numerical work as well as in pre- 
paring drawings. 
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Abstract 


The magnetic susceptibility of several olivines, «Fe,SiO,-(1—x)Mg,SiO, 
with varying x has been measured. The results showed that the olivines are 
paramagnetic and that their molecular magnetic susceptibility y,,., can be expressed 
empirically by ¥n9;=2x x10? °”"/.»., at the room temperature. Using the accepted 
values of the Bohr magneton number of Fet? ion, and assuming that the Fet? 
ions are embedded isolatedly in the crystal, the paramagnetic susceptibility of 
olivines at the room temperature can be calculated as ¥.,=2.2 x~2.6x x 107 
emt] 4,, provided that the other constituents such as Mg*? and SiO, are only 
diamagnetic. The comparison of the empirical and the calculated values seems 
to justify the above simple calculation. We expect such calculation will also valid 


for other rock-forming paramagnetic minerals. 


1. Introduction 


It is now an established fact that the magnetic properties of rocks relevant 
to geophysical phenomena are principally attributable to those of the ferromagnetic 
minerals contained in the rocks, such—as iron-titanium oxides and iron-sulphides. It 
does not seem to the present authors, therefore, the magnetic properties of so-called 
non-ferromagnetic minerals, such as silicates, which constitute over 95% of common 
rocks have ever attracted much attention of the investigators of rock-magnetism. In 
the present study, the magnetic susceptibility of some olivines, which are among the 
commonest rock-forming minerals, has been measured in order to clarify the para- 
magnetic nature of silicate minerals containing magnetic ions (Fe*t? in the case of 
olivines xFe,SiO,-(1—x)Mg.SiO,). 

Among the rock-forming minerals, ionic crystals which are composed of the ions 
with complete electronic configurations like those of inert gas atoms (!S)) and covalent 
crystals in which the valence electrons are antiparallel in spin are diamagnetic, while 
‘those minerals containing ions with incomplete shell configurations as the transition 
elements are presumed to exhibit paramagnetic properties.” Since olivines are the 
members of the solid solution between fayalite (Fe,SiO,) and forsterite (Mg,SiO,), in 
which the bonding between metallic ions and SiO, tetrahedrons is ionic and the bond- 
ing inside the tetrahedrons is covalent, the value of magnetic susceptibility is expected 
to vary with the ratio of forsterite to fayalite. mes 


-<* Contribution from Division of Geomagnetism and Geoelectricity, Geophysical Institute, Tokyo 
University, Series II, No. 66 ; 
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2. Samples Used and the Method of Measurement 


The samples examined in the present study, rendered by Prof. H. Kuno, Geologi- 
cal Institute, Tokyo University, are listed in Table I. The chemical composition of 
them was determined either by chemical analysis or measurement of the refractive 
index. 

After pulverized, each sample was fractioned by a Hallimond-type magnetic 
separator manufactured by Mitamura Co. Tokyo. Some results of the magnetic 
separation are illustrated in Fig. 1: the ordinate represents the amount of fractions 
yielded in traps A,B,C,D and E, which are in a decreasing order of the intensity of 
magnetism. These fractions will be called, in the following, the fraction A, fraction 
B, etc. Franction A, small in quantity, was always identified as magnetite through 
the measurement of the temperature variation of magnetization and by X-ray 
analysis. It was, moreover, observed that the fraction B could be sub-fractioned into 
the fraction A and the weaker fractions by further pulverization. Although these 
ferromagnetic impurities are small in quantity, they are the predominant factor in 
determining the apparent magnetic properties of the natural olivines, because the ferro- 
magnetism of magnetite is stronger than the paramagnetism of olivines by a factor 
larger than 10°*. 


Table I. Olivine samples examined and their magnetic susceptibility at room temperature. 


: Paramagnetic susce tibilit 
Locality Source Fa ee oe it in mol.| specific een 
y core (10-° emu/gr) |\(10-? emu/mol 
lizakamura, iron knebelite in 
Hukusima Pref. |pegmatite 98.4 10.3 20.8 
Huzi-Misima, (D)\olivine phenocrysts 
Sizuoka Pref, in basalt SBy 292? 2.7 4.3 
ak (E) ” 34~31 seul. 5.0 
Mt. Muie, (C)idunite in Bushveld 
Transvaal complex 33~30 2.3 3.7 
”» (D) » 34~29@) 2.5 4.0 
” (E) ” 32~29® Prat 4.0 
Isizumigahana, inclusion in anorthite 
Hatizy6-zima crystals in andesite 17.5@) 2.4 3.6 
Island pumice 
_|Itinomegata, olivine-rich inclusion 
Akita Pref. in basalt 5 10.6 1.8 2.6 
a ie (P)torsterite in marble 4@) 0.4 0.6 
” (E) 0 42) teh 1.6 


(C), (D), (E) mean the fractions separated by the magnetic separator. 

@K. Omori and S. Hasegawa, Sci. Rep. Téhoku Univ. ser. III, vol. IV, no. 1 (1951), 
@Private communication from S. Aramaki, Geological Institute, Tokyo University. 
@®Private communication from N. Issiki, Geological Survey, Japan. 

®C.S. Ross, M.D. Foster and A.T. Myers, Amer. Mineralogist, 39, 693-737 (1954). 


* Separation of minerals by magnetic separator such as Hallimond-type one is, thus, considered to 


be effected by rather the difference in the amount of ferromagnetic iopunities than the difference 
in the intensity of inherent paramagnetism. 
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Fig. 1 Examples of magnetic separation by a Hallimond-type magnetic 
separator: A....olivine in dunite, Transvaal: B... .olivine phenocrysts 
in basalt, Huzi-Misima: C....forsterite in marble, Miyako. 


Consequently, we studied the magnetism of pure olivines after getting rid of the 
magnetite impurities as perfect as possible. 

The magnetic measurements were all conducted by a sensitive quartz spring 
magnetic balance described elsewhere [1]. In spite of laborious removal of the magnetite 
impurities, even the weakest fraction of the olivines (fraction E) was found to contain 
a little amount of magnetite. In Fig. 2, curve a schematically shows the magnetization 
curve of a sample examined. The curve was interpreted as composed of the magnetiza- 
tion curve of paramagnetic olivine (curve 6) and that of ferromagnetic impurity (curve 
c). Therefore, the paramagnetic susceptibiliy of olivine was obtained from the slope 
of the straight line part AB. In the actual process, the effect due to the imperfect 
saturation of the ferromagnetic part was also taken into account. On certain samples, 
the temperature dependence of the susceptibility was also measured and it was ascer- 
tained that the initial curved part of the magnetization curve disappeared at temperatures 


higher than the Curie point of magnetite. 


Fig. 2 Schematic representation of magnetiza- 
tion curve of olivine: curve a@....curve 
actually measured: curve b....curve for 
purely paramagnetic olivine: curve C..+.- 
curve for ferromagnetic impurity. 


3. The Results of Measurement 


The molecular magnetic susceptibility at the room temperature of the samples 
measured by the above method is tabulated in Table I and plotted in Fig. 3 against 
the ratio of fayalite to forsterite. Some examples of magnetization curves measured 
are shown in Fig. 4, where the hollow circles represent the magnetization of the 
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fractions C,D and E of the olivine from dunite at Transvaal and the full circles are for 
the fraction E of the olivine from anorthite crystals in andesite pumice at Hatizyo-zima 
Island. In Fig. 4, it may be observed that the fractions C,D,E of the Transvaal sample 
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Fig. 3 Relation between the een susceptibility of olivine’ and the | 
; fayalite content at room temperature. 
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are nearly equal in para- 


as they differ in the amount of 
ferromagnetic impurity. This 
result seems to support the 
statement in the footnote in § 2. 

From the results in Fig. ; 
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the constituents. Therefore, the resultant magnetic susceptibility % can be expressed as, 


X=Aptxo; (2) 


where xp and x, representing the paramagnetic and diamagnetic susceptibilities. 
Although the contribution of diamagnetism in the present substance is not known, it 
may be presumed that we can negléct the second term in (2) to a good approximation. 

The effective Bohr magneton number pz of Fet? ion in crystal is known to be 
5.1~5.55 when calculated from the Curie’s law [2]. Using these values for the effective 
Bohr magneton number in the Curie’s law, the magnetic susceptibility per gramme ion 
of Fet’, x,.*? is enumerated as, 


pre Neer aoe 3.23 3.82 


3kT T iS 


(3) 


where N,,z, k and T are the Avogadro number, the Bohr magneton, the Boltzmann’s 
constant and the absolute temperature respectively. If, further, we assume that the 
magnetic moment of Fe*? ion is isolated from that of other Fe** ions in the olivine 
crystal, the molecular and the specific magnetic susceptibilities of fayalite (Fe,SiO,) at 
the room temperature become, 


Xmol fayalite — y Ea Ng ie 2.2 x 10-*~2.6 X 10-*emu/mol. 
% fayatttee = 11 x 10-°~13 x 10-* emu/gr. (4) 


As a result, the molecular magnetic susceptibility of any member of the solid solution 
xFe,SiO, « (1—x)Mg,SiO, at the room temperature can be expressed as, 


dno? ottvine = 2.24 x 10-2~2.64 x 10°? emu/mol. (5) 


Taking the experimental error introduced by the unavoidable presence of the 
ferromagnetic impurities into consideration, the agreement between the empirical 
formula (1) and the calculated one (5) seems to support the idea that the magnetic 
properties of olivine is. attributable to the paramagnetism of Fet*® ions isolatedly 
embedded in the crystal structure. If such agreement is to be expected to hold as a 
general rule, the paramagnetic susceptibility of various minerals will also be obtained 
by simple calculation of the following formula, 


2 
X mot = AEE Smbod = a x Si Pais (6) 


where the summation should be made for all the magnetic ions contained, and ”; and 
Pg, denote the number of the i-th magnetic ions in a molecule and the effective Bohr 
magneton number of the same ions. For instance, the molecular paramagnetic suscep- 
tibility of orthopyroxenes (xFeSiO3-(1 —*x) MgSiO;), namely the solid solution between 
orthoferrosilite (FeSiO;) and clinoenstatite (MgSiO;) can be expected to be the half of 


the value in (5). Further examination on various rock-forming minerals containing 


magnetic ions will be effective in confirming the present result as a general rule. 
Finally, the authors should like to express their sincere gratitude to Prot. i. 
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Kuno who kindly rendered them the useful samples and to Dr. S. Akimoto and Dr. S. 


lida for valuable discussions. 
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Disturbances in the Fl and FE Regions of the Ionosphere Associated 


with Geomagnetic Storms 


It was pointed out by Berkner and Seaton [1] that the critical frequencies (/,F1, 
fo) and the virtual heights ((h’F 1, h’E) of the F1 and E regions begin to vary at the 
same time when those in the F2 region do in the storm time. However, up to date 
the details of the variations in the Fl and E regions have been overlooked. One of 
the reasons may be that the variations in those regions at night cannot be known and 
the other that the deviations from the normal value are small compared with those in 
the F2 region. The purpose of this report is to show statistically the characteristics 
of the deviations in those regions. 

(i) In low sunspot or geomagnetic activities, the deviations of fF1l and f\k 
are small (especially the latter) and the amplitudes of them are often of the order of 
magnitude comparable with the error of the measurement. But with the increasing 
activity the amplitudes increase enough to admit as a disturbance. The disturbances 
in the Fl and E regions are classified in two types, the same as in the F2 distudance, 
positive and negative disturbances. The former is a disturbance in_which the critical 
frequency of ‘the region increases and the latter a disturbance in which it decrese. As 
shown in Figs. 2 and 3, the occurrence of the respective type of the disturbance depends 
on the season, latitude and geomagnetic or sunspot activities. The two types of the 
disturbance are shown in Fig. 1 for four stations. In that figure 4/F1 etc. represent 
deviative values from the respective normal values. Each variation is the mean of 
the variations on the first and second days for the moderate or severe disturbances. 
during the period shown in the figure. From Fig. 1 it is found that the daily variation - 
(only in daytime) of each type is similar to that in the F2 region. The deviation in 
the height seems to show an increase in the storm time. However, as Ah'E (sometimes 
Ah'F) is considered to be the magnitude comparable with the error in the height mea- 
surement and a correction of the retardation of the wave through the E region is 
not made for 4h’F1, the increments of the height in two regions are doubtful. 

(ii) In Fig. 2 the daily mean values (4/F1, AfE) of the 4fF1 and 4fE are 
plotted against the mean value (4f,F2) of the 4foF2 during the hours when the 4/,F1 
and 4fE are taken. From Fig, 2 the following points are recognized. 

(a) The positive correlation between the AF or AfpE and Af)F2 seems to be 
good. This correlation is better in high latitudes than in low and for 4foF1 than for> 

(b) The amplitude of 4/F1 and 4/E increases linearly in proportion to the’ 
amplitude of 4/oF2. 

(57 ) 
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Fig. 1 Mean daily variations of dfF1, dfE, 
Ah'F1 and dh’E for moderate and severe dis- (b) 


turbances during the period shown in the figure. Fig. 2 Relation between AfoF 1 and AfoE in 
Af F1 etc. represent the deviative values from the the disturbance occurred in four stations and 
respective normal values. The left-hand side shows the corresponding 4f,F2. dfF1 and AfoE 
the variation in the negative disturbance and represent the daily mean values of 4f)F1 and 
the right-hand side in the positive disturbance. AfpE and Af,F2 the mean value of AfoF2 
during the hours when f)F1 and /,E appear. 
The negative value represents the negative 
disturbance and the positive the positive dis- 
turbance. ° reperesents the value in summer, 
+ winter and * equinox. 


(c) In high latitudes foF1 and f,£ tends to decrease in the storm time (negative 
disturbance), in low latitudes (including the equator) they tend to increase (positive . 
disturbance) and in middle latitudes they increase generally in winter and decrease in 
summer and equinox. 

(iii) In Fig. 3 the values of 4fF and 4/,E for three stations are plotted against 
respectively YKp (sum of the three hourly Kp indices of a day) and against the pro- 


visional sunspot number at Zurich in the month when the disturbance occurs. The 
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Fig. 3 Relation between 4foFl and AfoE and =Kp (sum of three hourly Kp indices of a day) 
and the mean provisional sunspot number at Zurich in the month when the disturbance occurs. 
° represents fyF1 and* foE. 


ordinate represents 4/F1 and 4fk and the abscissa respectively SKp and sunspot 
number. From this figure the following points are recogniced. ; 


(a) AfoF 1 and 4foE increase in proportin to both Kp and sunspot number. 

(b) In middle latitudes the positive disturbance is easy to occur during the years 
of low sunspot or geomagnetic activities. 

(iv) From the above results it is found that the daily variation and the latitu- 
dinal and seasonal changes of the types of 4fpF1 and 4f,E in the storm time show the 
analogous feaures to those of AfoF2. These facts lead to a consideration that the 
disturbances in the Fl and E regions are caused by the same mechanism as that in 
the case of the F2 distubance, and it seems to be suitable to take the electron drift 
theory as the cause of the Fl and E disturbances as well as the F2 disturbance [2], [3]. 
Again, we are forced to admit the above theory the best from a fact that the disturbances 
in the F2, Fl and £ region begin. simultaneously, because the other theories, for 
example, the thermal expansion, recombination, are difficult to explain the above fact. 

From the view point of the drift theory we explain the statistical results shown 
above as follows (the same explanation is described in the case of the F2 disturbance 
[4], [5], so that the brief description is given here). In the storm time the electric field 
originated in the polar region spread over the world and the vertical drift of the elec- 
tron takes place due to this electric field and the geomagnetic. This velocity vz is 
dedueed from the geomagnetic disturbance-daily variation. On the other hand, the 
vertical drift with the velocity v¢ is present on a quiet day. The velocity is derived 
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from the Sg geomagnetic variation. On the disturbed day a vertical drift of the 
electron with the velocity (v¢+v,) occurs. The maximum electron density on a quiet 
day is, on:account of the presence of the vertical drift with the velocity v,, smaller 
than that in the static region. Therefore if the daily variation of vz has the same 
phase as that of vg, the reduction is promoted on the disturbed day. This is a negative 
disturbance. On the contrary, if the phase of the velocity v, is different by about 180° 
from that of v, the (vg+v.z) becomes samller than vg and the maximum electron density 
increases, because the region is analogous to the no-movement region. This is the 
positive disturbance. The former case is easy to occur in high latitudes and the latter 
in low latitudes (or in the equator), for the phase of v, differs depending on the latitude 
in spite of the invariability of the phase of v;. In middle latitudes the relation between 
the phases of vg and v, in summer and equinox is analogous to that in high latitudes 
and the relation in winter is analogous to that in low latitudes. Therefore the latitudi- 
nal and seasonal changes of the type of the disturbance appear as shown in Fig. 2. 
That the lower the region, the smaller the deviation is ascribed to the fact that 
the drift in lower region (much the same as in the upper region) becomes less effective 
due to the increase of the production and the decay of the electron. Martyn introduced 
the height gradient of the drift velocity in his rough estimation of the deviation in the 
F2 disturbance. Even if such a gradient exists, it is necessary that the gradients in the 
F2, Fl and E region are the same, inferring from the proportionality of the deviations 
as shown in Fig. 2. Actually, the necessary gradients of the drift velocity due to the 
electric field are probably difficult to occur in the Fl and E regions and also we do 
not need the gradient in the explanation of the F2 disturbance. Next, in the high 
sunspot activity the (Fl and f.E themselves (and also f)F2) are large compared with 
those in the low activity, so that the 4F1 and 4fE (4fF2) are larger, even if the 
deviation of the electron density is the same percentage. In the high geomagnetic 
activity (which is easy to occur in the high sunspot activity) the vz becomes larger 
compared with that in the low activity (that is, (vg+v,)~=vz). Then the variation 
(negative, or to the negative) is promoted since the drift with the veloety vg gives the 
effect of a reduction of the electron density in daytime, and somtimes removes the 
latitudinal and seasonal dependency. : 
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Meeting of the Society of Terrestrial Magnetism and Electricity : 
The 20th General Meeting was held at the Kyoto University on October 16-18, 
1956. 3 
Number of the Reports read af the Meeting: 
Rock Magnetism, 8; Atmospheric Electricity, 4; Radio Meteorology, 6; Geo- 
magnetism and Earth Current, 13; Ionosphere, 13; Cosmic Rays, 8. 
Tanakadate Prize was awarded for the following excellent worker : 
_ The 20th, Mr. T. Sato; 
Investigations on Disturbances in the Ionospheric F2 Region Associated with 
Geomagnetic Storms Based on the Ionization Drift Theory. 


The 21th General Meeting was held at the Tokyo University on May 10-12, 
1957. 
Number of the Reports read at the Meeting: 
Geomagnetism, 18; Ionosphere, 11; Night Airglow, 3; Cosmic Rays, 6; Rock © 
Magnetism, 9; Atomospheric Electricity and Radio Meteorology, 12; Special 
Report on Antarctic Research Expedition, 10. 
Tanakadate Prize was awarded for the oe excellent worker : 


The 21th, Mr. H. Ishikawa =e i> 
On waveforms of atmospherics of short distances and lightning: discharges 
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